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Abstract 

The possibility is considered for the formation in optical lattices of a heterogeneous 
state characterized by a spontaneous mesoscopic separation of the system into the spatial 
regions with different atomic densities. It is shown that such states can arise, if there are 
repulsive interactions between atoms in different lattice sites and the filling factor is less 
than one-half. 
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1 Introduction 


Optical lattices, loaded with cold atoms, are intensively studied, being the objects with rich 
properties that can be widely regulated (see, e.g., the review articles mi). The loaded atoms 
can interact with each other through short-range as well as long-range forces, such as dipolar 
forces [SHTJ. 

In addition to usual insulating and delocalized equilibrium states, atoms in optical lattices 
can form several quasi-equilibrium and metastable states. For example, in optical lattices, 
there can exist metastable states with repulsively bound bosonic pairs |8], metastable states 
characterized by microscopic phase separation in a mixture of two bosonic species |9], quasi¬ 
equilibrium mixture of localized and itinerant bosons HDl, and metastable states of atoms with 
dipolar interactions m- Double-well optical lattices can display the states with mesoscopic 
disorder, characterized by a heterophase mixture of mesoscopic regions with ordered and dis¬ 
ordered atomic imbalance [T2l[T3]. Incorporating into the system of cold atoms impurities [H] 
or imposing random external helds |15] can produce glassy lattice states ng similar to vitrihed 
solid states of metals na. 

In the present paper, we consider the possibility of forming in an optical lattice of a het¬ 
erophase state consisting of regions with different atomic densities. These regions have meso¬ 
scopic spatial sizes and are randomly distributed in space, where they are not hxed, but can 
appear and disappear in different places. In that sense, such a state is a dynamical heterophase 
mixture analogous to other heterophase states with mesoscopic phase separation, which occur 
in many condensed-matter systems [ 121111 ]. Each subregion of a competing phase is a kind of 
a droplet, or grain, of a denser phase inside a diluted phase. Such states are, of course, not 
absolutely equilibrium, but are quasi-equilibrium. 

The typical linear size of a dense droplet is dehned by the length /cor, at which atoms are 
strongly correlated and can coherently form a single phase. This length is mesoscopic, being 
between the mean interatomic distance a and the linear system size L, 

Cl Icor L . 


The droplet size is rather of nanoscale, not exceeding the critical radius, after which the germ 
would grow, provoking a phase transition in the whole system [20]. Nanoscale nuclei of a 
competing phase are not equilibrium and, strictly speaking, thermodynamic notions, such as 
surface tension or surface energy, may be not applicable [H]. The lifetime fcor of a correlated 
subregion, forming a droplet, is also mesoscopic, being between the local equilibration time tioc 
and the observation time tobsi 


tloc ^cor ^obs • 


Generally speaking, the sizes and lifetimes of the droplets are of multiscale nature, being inside 
mesoscopic intervals, for which l^or and Gor play the role of centers [ 22 ] • To some extent, the 
denser subregions remind the grains arising in the process of grain turbulence [23]. An opposite 
situation happens in the case of a solid with cracks and pores, where there are low-density 
regions inside a more dense solid [2^ . 

A snapshot of the heterophase two-density state is shown in Fig. 1, where the regions of 
higher density are randomly located inside a matrix of lower density. 

The consideration of a new thermodynamic state necessarily includes the analysis of its 
stability. Analyzing this, we show that there exist conditions, when the two-density state in 
an optical lattice is really stable. These conditions, briefly speaking, require the presence of 
intersite atomic interactions and a low hlling factor, smaller than one-half. 
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2 Heterophase two-density lattice state 


The first step for treating a two-phase system with random subregions is the averaging over 
heterophase conhgurations [T8l[T9]. Keeping in mind the standard form of the Hamiltonian, 
after averaging over conhgurations, we come to the effective Hamiltonian 

( 1 ) 


consisting of two terms 


H/y XCn 




Hdr) 


/i 


V’a(r) dr + 


tat 


V’I(r)t^l(r')<h(r - r')t/^„(r')'0„(r) drdr' (2) 


representing two different phases, whose atoms are described by the held operators 'ipa-, with 
a = 1,2. Here is an optical-lattice Hamiltonian and $(r) is a pair interaction potential. 

The Hamiltonian is renormalized by the geometric phase probabilities 

+ ^2 = K) , (3) 

where Va is the average volume occupied by the a - phase and V is the system volume. By this 
dehnition, the phase probability satishes the properties 


Wi + W2 = 1 , 0 < Wa < 1 ■ 

By assumption, the phases have diherent densities 

in which the number of atoms in an a - phase is 

Na = Wa I {'il^i{r)'ijja{r)) dr . 


( 4 ) 

( 5 ) 

( 6 ) 


Without the loss of generality, we may call the hrst phase more dense, so that 

Pi > P2 • (7) 

In that sense, the densities, distinguishing the phases play the role of the order parameters. 

The optical lattice prescribes the spatial periodicity of the lattice Hamiltonian with 

respect to the lattice vectors enumerated by the index j = 1, 2,..., running through all 
lattice sites. The held operators can be represented as expansions 

V'a(r) = (8) 

nj 

over the localized orbitals (pnj- The expansion takes into account that a j - site can be either 
occupied by an atom or free, depending on the value of the variable e“ = 0,1. 

We assume that each lattice site can host not more than one atom, which is expressed 
through the unipolarity condition 

= 1 , = 0 . (9) 


3 




Substituting expansion (8) into Hamiltonian (2) yields two types of terms, with respect 
to the site indices i and j. The terms, describing atomic interactions, dehne the effective 
time tosc of atomic oscillations in the vicinity of a given site. The other type of the terms is 
responsible for the hopping of atoms between the lattice sites, which can be characterized by a 
hopping time thop- The observation time has to be much longer than the hopping time, so that 
various phase conhgurations could be realized in the system, thus, justifying the averaging over 
these configurations. The relation between tosc and thop describes whether the system is in an 
insulating or delocalized state. When tosc is much shorter than thop, the atoms are well localized. 
This implies that the interaction terms are much larger than the hopping terms, responsible 
for atom hopping. In what follows, we assume that atoms are sufficiently well localized, so 
that the hopping terms are small, as compared to the interaction terms. Then the diagonal 
approximation can be employed corresponding to the following form of the matrix elements: 

{mi I Hl I nj) = SmJijEo , {mi,nj \ $ | m'i\n'j') = . (10) 

The constant term Eq can be incorporated into the chemical potential. In this way, Hamiltonain 
(2) reduces to the form 


Nl 


Nl 


Ha = -Wa^i e“ + 




J=1 




Using the canonical transformation 
1 




ej - 2 + (e“ - 0, l) , 


QZ _ a _ ^ 

^ ~ ^ 2 




we come to the pseudospin representation 




Nl 


2 Nl 


i=l 




where 


Nl 


Nl 


Then we resort to the mean-held approximation resulting in the Hamiltonian 


N 1 

Ha = [^2$ (1 - sl) - 4w„/i] + - [wl^ (1 + Sa) - 2Walj] ’ 

i=i 


(11) 

( 12 ) 

( 13 ) 

( 14 ) 

( 15 ) 


in which the notation 

2 

s„^2{y)„ = —(16) 

^ i=i 

is introduced. The average {Sj)a is taken with respect to Hamiltonain (15). Quantity (16) can 
be calculated either directly or by minimizing the thermodynamic grand potential 




a 



-2Tln2-TVln cosh 


a 


wl^{l + Sa) - 2WaH 

4T 


(17) 
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which gives 


Sq, 


tanh 


2vJaM - + Sq.) 

4T 


(18) 


Minimizing the grand potential (17) with respect to the phase probability, under the nor¬ 
malization condition (4), we set 


w = Wi , W 2 = 1 — w . 


( 19 ) 


The minimization yields the equation 


[2/iSa - w„<h(l 


dWa 

dw 


= 0 


from which we may express the chemical potential 

_ Wi(l -F Si)^ - W2{1 + S2)^ 

^ 2(si — S 2 ) 

An important quantity is the hlling factor 

1 ^=^ (TV = 7 Vi + 7V2) . 

tVL 

Representing the particle number (6) as 

Nt 

TVa = — Wa(l d- So) , 


4> 


we have the hlling factor 


V = 


y^wo(i+so) 


Inverting this with respect to the phase probability of the dense phase, we get 

2z/ — 1 — S 2 

Wi = - . 

-Sl ~ S 2 

It is convenient to introduce the dimensionless quantity 
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Xo = 


“ X] ’ 


i=i 


( 20 ) 


( 21 ) 


( 22 ) 


(23) 


( 24 ) 


(25) 


playing the role of a dimensionless order parameter. For the dense and diluted phases, we write 

(26) 


1 1 
X = Xi = -{l + Si), y = X2 = -{l + S2), 


respectively. The density of the a-phase reads as 

iVo Nr 




(27) 
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which shows why quantities (25) play the role of dimensionless order parameters. Due to 
inequality (7), we have the condition 

x>y , (28) 

distinguishing the phases with respect to their densities. 

The hlling factor (23) can be written as 

z/ = Wix + W 2 y ■ (29) 

Because of condition (28), the relation 

y <v < X (30) 


holds true. 

Using notation (26) reduces the chemical potential (20) to the form 

wix'^ - W2y‘^ , 

p =-<I> . 

x-y 

From Eq. (24), we find the phase probabilities 

ly — y X — V 


Wi = 


x-y 


W2 = 


x-y 


The sign 


$ 


a = — = sgn<F 


(31) 


(32) 


(33) 


of the effective interaction (14), for a while, is arbitrary. Measuring temperature in units of 
|$|, for the order-parameters (26), we get 


2x = 1 + a tanh 


2y = 1 + (T tanh 


wiy{wix - W2y) 
2T{x-y) 

W2x{wiX — W2y) 


2T{x-y) J ■ 

Thermodynamic quantities can be found from the free energy, for which we define the 
dimensionless quantity 


F = 


D -|- fiN a 


= I [wlx{l -x)+ wly{l - I/)] + Y ln[x(l - x)y{l “ ?/)] + y|j 


(35) 


3 Stability of heterophase two-density state 


First of all, we recall that to be stable a heterophase system has to satisfy the necessary 
heterophase stability condition 

which follows from the minimization of the grand potential [T81IT9] . This leads to the inequality 


/d^HX _ 


\ dw"^ 


= Nl^{x^ -f 2/^) > 0 , 


(37) 
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from which it is clear that the effective interaction (14) has to be effectively repulsive, so that 


$ > 0 


(38) 


An effectively attractive interaction does not allow for the formation of a stable heterophase 
system. 

Additionally, the system should be thermodynamically stable, implying that the specific 
heat 

'd^F\ 


Cv = -T 


dT^ J 


V 


and the isothermal compressibility 




) rj. 


be non-negative and finite 


0 < Cv < oo , 0 < kt < oo . 


(39) 


(40) 


(41) 


In what follows, we again measure temperature in units of <1>. Keeping in mind that we 
need to consider only repulsive atomic interactions, we have to solve the system of equations 
for the order parameter of the dense phase 


2x = 1 -I- tanh 
the order parameter of the rarefied phase 

2y = 1 + tanh 
and the probability of the dense phase 


wy[wx — (1 — w)y] 
2T{x-y) 

(1 — w)x[wx — (1 — w)y] 
2T{x-y) 


u-y 


w = 

x-y 

The corresponding solutions define the free energy 

1 1 
F = - \yx{l — x) -I- (1 — wYy{l — y)\ + -T ln[x(l — x)y{l 


(42) 


(43) 


(44) 


+ 


V — 


1 \ wx"^ — (1 — w)y‘^ 


(45) 


2/ X — y 

from where the specific heat and compressibility can be calculated. 

Solving Eqs. (42), (43), and (44), we are looking for the probability in the interval 0 < tc < 1 
and for the order parameters satisfying the inequalities 


Q < y < V < X <1 . 

Numerical investigation shows that the heterophase system can be stable only for small filling 
factors, 

0 <.< 1 . 
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For larger filling factors, compressibility (40) becomes negative, although specific heat (39) is 
always positive. 

Solutions to Eqs. (42) to (44) exist in the temperature interval [T„, T*]. These temperatures 
can be called the lower nucleation temperature and upper nucleation temperature. The lower 
nucleation temperature T„ is found numerically, being zero for v < 0.32. The upper nucleation 
temperature T* is defined by the conditions 

w(t:) = 0, x{t:) = i , y(T:) = i. , 

which yields 

^ (l-2zz)ln(i-l) • 

When z/ —)■ 1/2, then both and T* tend to infinity. Table 1 gives the values of and T* in 
the allowed interval of 0 < z/ < 1/2. 

Figure 2 presents the behavior of solutions as functions of temperature for different filling 
factors and Fig. 3, as functions of the filling factor for different temperatures. Specific heat 
(39) and compressibility (40) are positive. For illustration, Cy, as a function of temperature, 
is shown in Fig. 4. 


4 Comparison with pure single-density state 


The heterophase two-density state should be compared with the pure single-phase state, when 
w = 1. Then the grand potential is 


i <F(1 — s^) — \ jjii — Tin2 — Tin cosh 
Nl 8 2 


<h(l -1- s) - 2/ii 
4T 


where 


s = tanh 


2/ii — <h(l -|- s) 

4T 


The filling factor reads as 


N 1 

z/ = — = - (1 + s) . 
Nl 2^ ’ 


From Eqs. (47) and (48), we have 

2z/ = 1 -|- tanh 

which results in the chemical potential 


Pi — z/<h 
2T 


Pi = z/<F -I- T In 


1 - z/ 


For the dimensionless free energy, we get 


Ti = 


111 -|- pifV 1 2 


Nl\^\ 


= - + T[z/ In z/ -|- (1 — z/) ln(l — u) . 


( 46 ) 

( 47 ) 

( 48 ) 


( 49 ) 


(50) 



This shows that, for the pure phase, the specific heat is zero, Cy = 0. And for the compress¬ 
ibility, we find 

_ 1 ^ 

^ Z/[(TZ/(1 — v) +T] 

The latter is positive, provided that 

T > —au{l — v) . 

In the case of repulsive interactions, when a = 1, the pure phase can exist at all temperatures. 
But for attractive interactions, when a = —1, the system is stable only for sufficiently high 
temperatures, such that T > z/(l — u). 

Comparing the free energy (45) of the heterophase two-density state with the free energy 
(50) of the pure single-density state, we find that, in all those cases, when the heterophase state 
exists, F < Fi. This is shown in Fig. 5 for repulsive interactions. Therefore, in the temperature 
region Tn < T <T*, the heterophase state is stable, while the pure state is metastable. 

5 Conclusion 

We have considered the possibility of the formation in optical lattices of a heterogeneous state 
characterized by a spontaneous mesoscopic separation of the system into the spatial regions 
with two different atomic densities, one being more dense than the other. We show that such 
states can really occur, provided that atomic interactions between atoms in different lattice 
sites are repulsive and the filling factor is less than one-half. The heterophase state is stable in 
the temperature region between the lower, T„, and upper, T*, nucleation temperatures. 
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Figure Captions 


Figure 1. Snapshot of a heterophase two-density lattice system. Regions of higher density 
Pi are randomly immersed into the matrix of lower density p 2 , with pi > p 2 - 

Figure 2. Solutions as functions of dimensionless temperature T for different hlling factors: 
(a) order parameters x (solid line) and y (dashed line) for u = 0.1 (line 1), u = 0.2 (line 2), 
u = 0.3 (line 3), and u = 0.32874 (line 4); (b) dense-phase probability w for the same hlling 
factors and enumeration as above; (c) order parameters x (solid line) and y (dashed line) for 
z/ = 0.33 (line 1), z/ = 0.4 (line 2), and z/ = 0.45 (line 3); (d) dense-phase probability w for the 
same hlling factors and enumeration as in (c). 


Figure 3. Order parameters x (solid line) and y (dashed line) and the dense-phase prob¬ 
ability w (dashed-dotted line) as functions of hlling factor v for diherent temperatures: (a) 
T = 0.01; (b) T = 0.5. 


Figure 4. Specihc heat as function of temperature for diherent hlling factors: (a) z/ = 0.1 
(line 1), z/ = 0.2 (line 2), z/ = 0.3 (line 3), and u = 0.32874 (line 4); (b) z/ = 0.33 (line 1), 
u = 0.4 (line 2), and z/ = 0.45 (line 3). 


Figure 5. Free energy F (solid line) of the heterophase state, compared to the free energy 
Fi (dashed line) of the pure state, for varying temperature and diherent hlling factors: (a) 
z/ = 0.3; (b) z/ = 0.45. 
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Table Caption 


Table 1. Lower and upper temperatures defining the existence interval of the optical lattice 
with heterogeneous densities. 


Table 1 


V 

T 

n 

0.1 

0 

0.05689 

0.2 

0 

0.24045 

0.3 

0 

0.88517 

0.32874 

0.01 

1.34442 

0.33 

0.01253 

1.37053 

0.4 

0.520515 

4.93261 

0.45 

3.9583 

22.4248 
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Figure 1: Snapshot of a heterophase two-density lattice system. Regions of higher density pi 
are randomly immersed into the matrix of lower density p 2 , with pi > p 2 . 
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Figure 2: Solutions as functions of dimensionless temperature T for different filling factors: (a) 
order parameters x (solid line) and y (dashed line) for v = 0.1 (line 1), v = 0.2 (line 2), v = 0.3 
(line 3), and v = 0.32874 (line 4); (b) dense-phase probability w for the same hlling factors and 
enumeration as above; (c) order parameters x (solid line) and y (dashed line) for v = 0.33 (line 
1), z/ = 0.4 (line 2), and v = 0.45 (line 3); (d) dense-phase probability w for the same hlling 
factors and enumeration as in (c). 
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Figure 3: Order parameters x (solid line) and y (dashed line) and the dense-phase probability 
w (dashed-dotted line) as functions of hlling factor v for different temperatures: (a) T = 0.01; 
(b) T = 0.5. 
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Figure 4: Specific heat as function of temperature for different filling factors: (a) i/ = 0.1 (line 
1), V = 0.2 (line 2), v = 0.3 (line 3), and v = 0.32874 (line 4); (b) v = 0.33 (line 1), v = 0.4 
(line 2), and v = 0.45 (line 3). 
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Figure 5: Free energy F (solid line) of the heterophase state, compared to the free energy Fi 
(dashed line) of the pure state, for varying temperature and different Filing factors: (a) v = 0.3; 
(b) z/ = 0.45. 
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